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ABSTRACT 


In the present paper the fitting and proper regression coefficients have been made of 
117 1 0'xl O' blocks with observed gravity data and corresponding elevation in the Taiwan 
Island To compare five different predicted models, and the proper one for the mean 
gravity anomalies were determined. The predicted gravity anomalies of the non-observed 
gravity blocks were decided when the coefficients obtained through the model with the 
weighted mean method. It was suggested that the mean gravity anomalies of IO’xIO' 
blocks should be made when comprehensive the observed and predicted data. 

1. INTRODUCTION 


The purposes of the paper is to understand the relationship of gravity anomaliy with 
the topography and the area, doing further research, comparing the good and bad points 
from the most common used mean computing models. Analysising which model is more 
suitable for Taiwan Island and meatime computing the accuracy of mean gravity anomaly. 
This paper is mainly stressed on Taiwan area, surrounded by the longitude from 120° E 
to 122° E, the latitude from 21.5° N to 25.5° N. 


2. COMPUTATION MODELS 


There are five formulas of the mean gravity anomaly will be discussed [Uotila, 
1967a&b ; Kassim, 1980;Sunkel, 1983 ]: 


(1) Ag = a + bK , 

(2) 5g = 1/n XAg; _ 

(3) Ag = 1/n I [Agj - b (hj-H)] t 

(4) Ag = I (Agj /sj 3 5 )/X(1/Sj 3 5 ) , 

(5) Tg = {X [(Agj -bhj)/Si 3 ' 5 ]/I(1/Sj 3 ' 5 )} + bh , 
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which Ag is the mean gravity anomaly of square block; Agj is the ith gravity anomaly in 
square block; h is the mean height of square block; hj is the ith height in square block; a, 
b are the regression coefficients; Sj is the distance between number the ith point and the 
centre of the square block; 3.5 is the weighted exponent. 


In these five computing models, the needed parameters as Agj, hj, h, Sj, a, b. And Agj, 
hj are obtained from the observation values, Sj is the distance betwen the observation 

point and the centre point in block, h is obtained from digital terrain model. Therefore, 
the regression coefficients a, b are obtained from the first order of the stochastic 
functions of the gravity anomaly and height. In this paper, they come from two mam 
resources: (1) using more than two points data in every individual block to calculate its 
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a, b value in block, (2) using all observed gravity materials, setting different gr'ou'ps 
according to the height of topography then calculate the a, b value in three different 
topography in Taiwan area (Table 1). 

Table 1 . The three groups of regression coefficient value of heights 


HEIGHT(m) 

a (mGal) 

b(mGal/m) 

h < 100 

100 < h <1000 j 

h > 1000 

- 1.8585 
13.7330 
-63.8362 

0.1091 

0.0281 

0.1195 


Therefore, as long as using mean height in each block, the selected a, b coefficient can be 
determined, a, b are calculated with data from all over Taiwan, so the stability is very 
strong. When the height is between 100 metre and 1,000 metre, the value of b is 
regarding 0.0281 to 0.119, somewhate different in theory but it is due to the effect of 
topography and area. 


3. THE RESULT AND ANALYSIS 

Having 117 blocks in observed gravity data, it makes difficult to show them all. We 
picked 12 block results at random to discuss as a base line. We select different 
topography - plain, hill, mountain besides chosing each block in different location from 
the east to west and the south to the north. In keeping with comparing and analysing the 
effect from five different computing models, the results of above selected 12 blocks is 
showed in Table 2, and among the five numbered gravity are corresponding with the 
sequence of computing data. 

As following, we point out and explain some results, we compared and analysed: (1) 
Clearly, the computed result from five models may be devided into three categories: 5g-| 
as a independent, Ag 3 , Ag4 as category, Ag 3 , Agg is another one. Ag2 - ^94 as deficit 
result of height and gravity anomaly which the relationship between them has not been 
considered. And £g 3 , Ag 5 as a group result from the effect of height was considered. 
Then why the result is differentin Ag-| with the effect of height? The reason in Ag-, = 
a+bh within, we set a as fixed value. In the contrary, by using a value for section area to 
compute 5g-| , the result is unified in mathematical meaning and Ag 3 (2) The different 
between the group ( Ag 3 , Ag4) and Ag 3 Ag 5) is explained by location of height. The 

difference in these two depend on if there is affecting existance from height to gravity 
anomaly. Thus inter-relation may very well became greater by increasing the height. In 
the contrast, the deficit is greater comparing mountain with hill, the hill somewhate 
larger than the plain. Therefore the relationship between height and gravity anomaly 
have to be stressed on particularly the mountains areas. Providing standardizesd data in 
mean gravity anomaly in any topography, we must disregard height in predicted model. 
Then the group of Ag 3 and Sg 4 is abandoned. (3) Finally, in group Sg 3 and Sg 5 , the 

reason, the deficit existed is if the locations of points can represent the mean position as a 
whole in any blocks. i^g 3 is computed from any point position, 5g5 is using the position 
of the centre point in the block to compute mean value. The closer the points to the 
centre are, the closer result of 5g 3 and ^5 are. Again if the points in block are 
situated uniformly, the Ag 3 value will be the most accurate. In 5g 3 , it is no way to 

prove if the mean gravity in centre position can represent mean value according to the 
result obtained from Table 2, there are a small discrepancy and Ag 3 computing model is 
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simpler and the more idealistic positions are, the more reasonable result are. 

Resent paper decided to use Sg 3 as computing method for mean gravtty anomaly 


model. 


Table 2. The results of five computation models 


NO^XMEAN 

OF \ GRAVITY 
BLOCK \ ANOMALY 

Agl 

Ag2 

&g3 

Ag4 

Ag5 

179 

183 

63 

-1.86 

6.87 

3.60 

-31.30 

-25.36 

-14.39 

-58.66 
-17.81 
- 9.79 

-31.3 

-24.67 

-14.03 

-58.66 

-16.76 

-9.43 

260 

165 

185 

3.60 

36.78 

24.97 

2.69 

50.64 

-37.84 

- 4.55 
69.65 
-29.51 

13.04 

38.83 

-41.60 

4.77 

58.17 

-34.08 

53 

250 

177 

29.19 

28.49 

95.10 

65.98 

9.76 

39.10 

77.14 

20.25 

166.73 

66.6 

17.89 

34.68 

77.23 

25.56 

159.05 

162 

113 

237 

121.39 

68.81 

82.55 

-17.88 

19.47 

38.59 

122.08 

96.11 

122.84 

-17.88 

18.12 

36.02 

1 22.08 
95.41 
117.97 


4. THE ACCURACY OF THE MEAN GRAVITY ANOMALY 

Because of being unable to observe the exact mean gravity anomaly in each block, we 
Jfe no way to discuss the experimental accuracy of the block, we can use the law of 
propagation Of errors to eshmate the accuracy of the block. The present paper deeded ,0 
use the mean gravity anomaly computing model, that 


Ag = 1/n I [Agj - b (hj-h)] i 

we can predict mean gravity anomaly in block in Taiwan Island with observed data (Table 
3). 


Table 3 Categorlized by topography to determine the 
accuracy of the mean anomaly in the block 


' ---~-___TOPOGRAPHY 1 

accuraoy~~~ — — — _ 

PLAIN 

HILL 

mountain 

a Ag (mGal) 

4.8 

3.5 

~~ 16.1 
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5. CONCLUSION 


Being the 2/3 of Taiwan Island is mountain, situated on earthquake zone, the gravity 
anomaly is greatly different by year. At present, Taiwan Island we have just completed 
the levelling stations to promote gravity observations. Along with doing dense in 
mountain area, we strongly believer we can obtain more accurate mean gravity anomaly 
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